A room temperature 19-channel magnetic field mapping device for cardiac signals 
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We present a multichannel cardiac magnetic field imaging system built in Fribourg from optical 
double-resonance Cs vapor magnetometers. It consists of 25 individual sensors designed to record 
magnetic field maps of the beating human heart by simultaneous measurements on a grid of 19 points 
over the chest. The system is operated as an array of second order gradiometers using sophisticated 
digitally controlled feedback loops. 



Signals derived from the electrophysiological processes 
in the human heart are important diagnostic tools for as- 
sessing cardiac disorder. For more than a century, mea- 
surements of body surface potentials (electrocardiogra- 
phy, EGG) generated by electrical currents within the 
myocardium have been a standard medical examination 
technique. Electrophysiological currents also generate a 
magnetic field, whose recording represents an alterna- 
tive diagnostic tool known as magnetocardiography 
(MCG). Three decades after the first recording of an 
MCG signal 0, Schneider et al.^ in 1990 turned the 
MCG method into an imaging technique (cardiac mag- 
netic field mapping, cMFM) by using a sensor array cov- 
ering the whole chest. The magnetic signals (MGG) have 
a distinct advantage over the electrical signals (EGG) 
when considering source localization: the inverse prob- 
lem of reconstructing the field sources from field maps 
is easier to solve in the magnetic case since the electri- 
cal conductivity of body tissue affects the transmission 
of magnetic signals to a much lesser degree than their 
electric counterpart. Modern cMFM systems [3| use su- 
perconducting quantum interference devices (SQUlDs) 
for recording the very weak magnetic field (peak ampli- 
tude ~ 100 pT) generated by the heart. In the past years 
we have developed laser optical pumping magnetometers 
(LOPMs) [H for biomagnetic applications @ . The room 
temperature operation of LOPMs makes such magne- 
tometers promising alternatives to SQUIDs by avoiding 
the expensive cryogenic cooling and the associated com- 
plex logistics in view of clinical applications. Since each 
LOPM sensor requires only a few jiW of light power, a 
single laser can operate many dozens of individual sensor 
heads, a feature which has until now been unexploited. 
In this Letter we describe a cMFM system operating in 
Fribourg consisting of 25 individual LOPM sensors de- 
signed to record magnetic field maps of the beating hu- 
man heart by simultaneous measurements on a grid of 
19 points over the chest. 

The LOPM sensors are based on optical-rf double res- 
onance in the so-called configuration Q . The sensor 
medium is a room temperature Gs vapor contained in a 



spherical paraffin-coated Q 30 mm diameter glass bulb. 
Gircularly polarized light from a diode laser, frequency 
locked to the 4-^3 component of the cesium Di transi- 
tion, creates a macroscopic magnetization in the vapor by 
optical pumping. The magnetization precesses around a 
static magnetic field Bq = Bq z — nominally oriented at 
7r/4 with respect to the light propagation — with the Lar- 
mor frequency ujl = 7 |i3o|, where 7 « 27r x 3.5 Hz/nT. 
The precession is coherently driven by a weak magnetic 
field Bi ^ Bq oscillating at the frequency Wrf ■ As a con- 
sequence, the light transmitted through the cell acquires 
an amplitude modulation at the frequency cj^f • The mod- 
ulation amplitude is resonantly enhanced at coj-f = lul and 
the phase A(/3 between the drive and the response obeys 
(up to a constant phase offset) 



= arctan 



(1) 
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where the (HWHM) linewidth F is typically « 27r x 6 Hz 
for our OPMs. Near resonance, depends linearly 
on Bq and on Wrf, so that feedback control of either Bq 
or Wff can be used to keep the system on resonance. 

The design of the 19-channel system was targeted 
at demonstrating the simultaneous operation of many 
LOPM sensors in close packing. Each sensor consists of 
a compact self-contained module including optical com- 
ponents to coUimate and polarize the light, a Gs vapor 
cell, and a photodiode to detect the light transmitted 
through the cell. The components arc mounted between 
two printed circuit boards (PGB) which carry coils for 
generating Bi. Each sensor module is placed in an ad- 
ditional pair of coils allowing the application of a small 
field, Sfb, parallel to the field Bq. These feedback coils 
are realized as printed circuits on larger PGBs that serve 
as mechanical support structures positioning the sensors 
in horizontal planes (Fig. [T]). An optimal sensor pack- 
ing is achieved by arranging the modules on a hexagonal 
grid with a 50 mm spacing. A custom-made hologram 
splits the beam from the frequency-stabilized laser into 
25 beams of equal intensity which are carried by mul- 
timode optical fibers to the individual sensor modules. 
The complete sensor system consists of 25 modules dis- 
tributed in three planes that are vertically offset by 10 cm 
(Fig. [T]). The lower plane, closest to the subject's chest, 
is equipped with 19 primary mapping modules. The mid- 
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Figure 1: Perspective drawing of the 25-sensor 19-channel 
cMFM system installed in our Al shielding room in Fribourg. 
The bottom plane holds the 19 measurement channels, while 
the mid plane and top plane hold three sensors each, used for 
field and gradient fluctuation suppression. 

die and upper planes each carry three sensors providing 
reference signals for stabilizing Bq and its linear gradi- 
ent, respectively. The sensor system is mounted inside 
of a two layer Al shielding room that reduces ambient 
magnetic field fluctuation which arc dominated by 50Hz 
(power line) and 16 2/3 Hz (railway) interference. The 
person to be diagnosed is placed underneath the sensor 
system on a nonmagnetic bed. Three pairs of large coils 
surrounding the shielding room are used to create a verti- 
cal field -Bo ~ 9/xT (corresponding touL = 27rx 31.5 kHz) 
at the sensor location. Transverse dBo/dx, dBo/dy gra- 
dients are compensated using additional large coils and 
constant current supplies. 

The control electronics is implemented in a digital 
FPGA board which digitizes, via ADCs, the photocur- 
rent from each module after transimpedance amplifica- 
tion, and which generates, via DACs, the currents driv- 
ing the rf coils and the feedback currents for the sensors 
in the lower and upper planes. For each sensor the phase 
difference Aip is determined using a numerical lock-in de- 
tection algorithm. The near resonant linear dependence 
of the phase signals form suitable error signals that are 
used, after proportional-integral-differential (PID) am- 
plification, as feedback signals in the various servo-loops 
discussed below. The middle plane carries three sensors 
that serve as reference sensors which track changes of 
the DC offset field Bq (applied field plus external field 
changes). In each of the three middle plane sensors, A(p 



is actively locked to zero by servo adjustments of the 
corresponding three drive frequencies uj,-[ generated by 
a numerically controlled oscillator. The computed aver- 
age of those three individual frequencies or, alternatively, 
the frequency of a single (centered) sensor is then used 
as common drive frequency LUcom for the sensors in the 
top and bottom planes. In both ways the drive frequency 
tracks fluctuations and drifts of Bq. 

When using a common drive frequency cjcom, the sig- 
nals Aip of the sensors in the top and bottom planes are 
proportional to the difference between the field at their 
locations and the field in the middle plane. These sensors 
thus form first-order hardware gradiometers. In particu- 
lar, the signal Atp of the top sensors measures the linear 
vertical magnetic field gradient, while being almost in- 
sensitive to the cardiac field. The phase differences Acp 
of the upper sensors are locked to zero via servo loops 
injecting current into each of their respective feedback 
coils. In this way the local field at each upper sensor 
location follows the field in the middle plane. When si- 
multaneously injecting currents of the same magnitude, 
but of opposite polarity, into the bottom layer feedback 
coils, the linear vertical field gradient along the vertical 
extension of the sensor system, and its fluctuations are 
actively compensated. Since the number of upper refer- 
ence sensors is less than the number of imaging sensors 
in the bottom plane, only the average of the upper cor- 
rection currents is used as feedback for the lower plane. 

When the average field and the linear vertical gradi- 
ent fluctuations are suppressed by the feedback loops de- 
scribed above, the remaining phase differences Aip of the 
lower plane sensors arc proportional to the second or- 
der gradient of Bq. However, those phase differences are 
not directly measured: by using yet another set of servo- 
loops driving the lower plane feedback coils, the 19 phase 
differences Aip are removed, forcing all sensors to be in 
resonance with c^Jcom- The 19 time-dependent feedback 
currents of the latter loop thus consist of the biomag- 
nctic signals of interest and fluctuations of Bq^s second 
order gradient. These signals are recorded at a sampling 
rate of 5 kHz. This magnetic feedback scheme relying 
on three distinct sets of servo-loops is of general inter- 
est for any gradiometric arrangement of multiple, closely 
packed magnetometers. It may be particularly useful for 
coherent population trapping magnetometers (CPTM), 
which have been used for biomagnctic field detection Q 
but have not yet been reported to be scalable to many 
channels. Close-packed CPTMs can be operated in an 
all-optical mode by using one optical modulator per chan- 
nel, however, magnetic feedback is probably a more cost 
effective operation method for a multichannel CPTM im- 
plementation. The second order gradiometer scheme re- 
duces field fluctuations by a factor of >1000 to a level 
of 900 fT/\/Hz at 0.1 Hz. The noise of the gradiometer 
signal drops towards larger frequencies and reaches 300 
fT/Vlfe at 100 Hz. 

Since the sensors are closely packed, the i?tb field of a 
given feedback coil induces crosstalk with its neighbours. 
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Figure 2: Time series showing (averaged) MCG traces of a 
healthy subject for each of the 19 sensor positions. Each 
trace is 1 s long and all graphs have the same (280 pT) verti- 
cal extension. For this measurement only a single (centered) 
reference sensor was used in the middle plane. 
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Figure 3: (Color online) Magnetic field map for the out-of- 
chest component of the cardiac field. Crosses show sensor 
positions, and the superimposed trace (bottom center trace 
from Fig. [2]) identifies, via the black dot, the time for which 
the map was evaluated. Left (right): maps at the R-peak 
(T-wave). 



although the effect is small enough (<10%) that the sys- 
tem operates without online corrections. The measured 
crosstalk dBi / dij (field in sensor i induced by current in 
feedback coil j) can be used for offline data correction]^. 

Figure [2] shows an example of 19 simultaneously 
recorded MCG traces obtained by averaging signals from 
66 heartbeats. The data represents the second order gra- 
dient of the out-of-chest component of the cardiac field. 
A 3-lead ECG signal served as timing reference for align- 
ing the individual heart beats when averaging. The total 
recording time is less than two minutes, a huge improve- 
ment compared to the two hours required when stepping 
a single sensor over the chest Q. The signals are (low- 
pass) filtered offline by multiplying their Fourier spectra 
with a smooth function passing from 1 to in the band 
of 40 to 50 Hz. Systematic implementation of crosstalk 
correction is a work in progress and thus has not been 
applied to the shown data. The residual oscillations are 
due to railway interference (16 2/3 Hz). In Fig. [31 mag- 
netic field intensity maps for two given times (R-peak, 
T-wave) of the cardiac cycle are shown. The smooth 
pattern is obtained by a numerical interpolation between 
the discrete measurement points (shown as crosses). 

We have reported the operation of a 19 channel second- 
order gradiometer system based on atomic magnetome- 
ters for MCG measurements, capable of recording a mag- 
netic field map in about two minutes in a weakly shielded 
environment. The 19 + 3 + 3 system is now being used 
at the University of Fribourg for an initial series of mea- 
surements on both volunteers and cardiac patients. A 
57 + 13 + 13 channel system is currently being built at 
the University Hospital in Jena (Germany) . Tests of that 
system in a clinical environment should start early in 
2010. 
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